Gravitational waves from inspiraling compact binaries: The quadrupole-moment term 
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A rotating star's oblateness creates a deformation in the gravitational field outside the star, which 
is measured by the quadrupole-moment tensor. We consider the effect of the quadrupole moment 
on the orbital motion and rate of inspiral of a compact binary system, composed of neutron stars 
and/or black holes. We find that in the case of circular orbits, the quadrupole-monopole interaction 
affects the relation between orbital radius and angular velocity, and also the rate of inspiral, by a 
quantity of order (n/c) 4 , where v is the orbital velocity and c the speed of light. 
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I. INTRODUCTION AND SUMMARY 



Inspiraling compact binaries, composed of neutron 
stars and/or black holes, are the most promising source 
of gravitational waves for kilometer-scale interferometric 
detectors such as the American LIGO Q and the French- 
Italian VIRGO j| . The current construction of these de- 
tectors has motivated a lot of recent work on strategies 
to detect and measure such waves . In the course of 

this work, the need for extremely accurate predictions of 
the expected signal has repeatedly been demonstrated, as 
reliable model signals will be required for analyzing the 
data with the well-known technique of matched filtering 

& 

A large effort is currently underway to calculate in- 
spiraling binary waveforms to high order in a post- 
Newtonian approximation to the exact laws of general 
relativity p3|-|l5|. This method is based on the assump- 
tion that the orbital motion is sufficiently slow that the 
waves can be expressed as a power series in v, the orbital 
velocity. Thus far, the waveforms have been calculated 
to order v 5 beyond the leading-order expressions, under 
the assumption that the spinning motion of the binary 
companions can be ignored |14|. Rotational corrections 
to the waveforms have been computed only up to order 

Of particular importance for matched filtering is the 
phasing of the waves [fj"7| , which is determined by the rate 
of increase of the gravitational-wave frequency / (equal 
to twice the orbital frequency) in response to the system's 
loss of energy and angular momentum to gravitational 
radiation. Assuming, as is usual, that the orbital motion 
is circular, this is given by || 



0{v 5 ) 



(1) 



where we use units such that c = G = 1. We have 
introduced a number of symbols. Let mi and 771,2 denote 
the masses of the two companions, and let M — mi +7772 
be the total mass and /x = m\m%lM the reduced mass. 
Then we define the chirp mass Ai and the mass ratio 77 
as 



M = rf /5 M, 77 = fi/M. 



(2) 



The orbital velocity v is defined in terms of the 
gravitational-wave frequency / by 



v = (nMf) 1/3 . 



(3) 



Finally, the quantities f3 and a have to do with the spin 
angular momenta of the companions, denoted Si and «2- 
The "spin-orbit" parameter is given by 



^ = 1^ E [H3(m A /A/) 2 + 7577] L ■ X A, 



(4) 



where L is the direction of the orbital angular momen- 
tum, and xa — ^a/it^a 2 is the dimensionless spin of 
companion A. (The index A runs over the values 1 and 
2. Here and throughout the paper, vectors and tensors 
are defined in three-dimensional flat space, and a hat in- 
dicates that the vector has a unit norm.) On the other 
hand, the "spin-spin" contribution to a is |l6| 



a ss = ^ [-247(xi ■ X2) + 721(i • X i)(L ■ X2)] ■ (5) 
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Our purpose in this paper is to derive an additional 
contribution to <r, namely c qm , which is due to the 
quadrupole moments of the companions. The physical 
picture is the following. The spinning motion of com- 
panion A creates a distortion in its mass distribution 
which, in turn, creates a distortion in the gravitational 
field outside the star, measured by Q A b , the (trecefree) 



1 



quadrupole-moment tensor. The quadrupole term in the 
gravitational potential affects the orbital motion of the 
companions, and it affects also the emission of gravita- 
tional waves. This is the effect that we consider in this 
paper. It is important to understand that we are not 
considering the gravitational waves emitted by the time 
variations of the quadrupole moments Q\ h . Although 
such an effect exists (because of spin precession), it is 
much weaker (by an estimated factor of order v 16 !) than 
the one considered here. 

Assuming that the spinning body A is axially symmet- 
ric about the direction of §a, the quadrupole-moment 
tensor can be expressed as 



I A 



a 



(6) 



where Qa is the quadrupole-moment scalar. In New- 
tonian theory fl8[ ], this is given in terms of the mass 
density p by Qa — J A p( x ) \ x \ 2 ^(S ■ x)d 3 x, where 
P2(x) = \{ix 2 — 1). In general relativity, Qa is defined 
in a coordinate-invariant manner in terms of the falloff 
behavior of the metric outside the star (see Ref. Jl?|| and 
references therein). The general relativistic definition re- 
duces to the Newtonian one when the gravitational field 
is weak everywhere inside the star. 

In the following sections of this paper we calculate that 
the quadrupole-moment contribution to a is 



'qm 



2 ^ m A M 2 



3(L • § A ) 2 - 1 



(7) 



reflects the oblateness of the black hole.) Quadrupole mo- 
ments of realistic rotating neutron stars have been com- 
puted in Ref. p9|, where it is shown that for a neutron 
star of 1.4 Mq, 



q, 



-axA 2 m A 3 - 



(8) 



Here, the parameter a ranges from approximately 4 to 
8 depending on the equation of state for neutron-star 
matter. Stiffer equations of state give larger values of a, 
and a = 1 for a rotating black hole. Inserting Eq. (||) into 
Eq. ^ shows that whatever the nature of the compact 
object (black hole or neutron star), o~ ss and <7 qm are of the 
same order of magnitude. For example, for two 1.4M Q 
neutron stars with a — 5.0, we find cr ss £ (—1.04,1.04), 
while <7 qm <E (—2.64,5.28). We have used the fact that 
for such neutron stars, \ < 0-65 |]l9| . 

The rest of the paper is organized as follows. In Sec. II 
we integrate the Newtonian equations of motion for the 
centers of mass of the two companions, focusing on the 
motion over a time scale comparable with the orbital 
period. In Sec. Ill we consider the precessional motions 
of the orbital and spin angular momenta, which occur 
over a time scale much longer than the orbital period. 
Finally, in Sec. IV we calculate the contribution to the 
rate of change of the gravitational-wave frequency due to 
the quadrupole-monopole interaction. 



II. ORBITAL MOTION 



where Qa is the general-relativistic quadrupole-moment 
scalar. This is the main result of this paper. To the best 
of the author's knowledge, this contribution was never 
presented in the literature before, except in the specific 
context of a particle moving around a massive Kerr black 
hole pop . We note, however, that similar calculations 
were carried out independently, but not published, by 
Kidder glj. 

The quadrupole-monopole interaction responsible for 
CT qm is a Newtonian effect which formally takes the ap- 
pearance of a second post-Newtonian correction in df /dt. 
Indeed, apart from the facts that Qa is defined in a 
general-relativistic manner and the Einstein quadrupole 
formula is used to compute the radiation, the calcula- 
tion leading to Eq. (0) involves only Newtonian theory. 
In other words, although the gravitational field is not 
assumed to be weak inside the compact objects, the rela- 
tive separation of the companions is assumed to be suffi- 
ciently large that the mutual gravitational potential will 
be small. Under such conditions, Newtonian theory can 
be used to describe the orbital motion, provided that the 
quadrupole moments are defined in a suitable, strong- 
field manner. 

It is well known that the quadrupole moment of a ro- 
tating black hole is given by Qa — — XA 2 ^riA 3 , where 
Xa = \xa \ is the dimensionless spin (p2|. (The minus sign 



The Newtonian equations of motion for the centers 
of mass of two bodies with masses mi and and 
quadrupole moments Q\ and Q2 can be derived from an 
effective one-body Lagrangian, C — \pLX 2 — V(x), where 
fi is the reduced mass, x = X2 — x\ the relative sepa- 
ration of the centers, and x the relative velocity. The 
gravitational potential is [E3| 



V(x) 



fiM 3 



(miQf + rniQf ) 



(9) 



where M is the total mass, r = \x\, n = x = x/r, and 
Q A b is constructed from Qa as in Eq. (||). 

For what follows it is useful to define the dimensionless 
quantities 



PA 



Qy 



tuaM 2 

and to introduce the angles and Pa such that 
§a = (sinoi^ cos Pa, sina^ sin/3,4, cosoia)- 



(10) 



(11) 



Although these angles vary because of spin precession 
(see Sec. Ill), this variation occurs on a time scale much 
longer than the orbital period. We shall therefore take 
them to be constant for the purpose of calculating the 
motion over a time scale comparable with the orbital 
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period. We also introduce the spherical coordinates 
{r, 0, (j)}, such that x = rn and 



n — (sin 9 cos 4>, sin sin (j>, cos t 



(12) 



Substituting these relations, together with Eq. (^|), into 
Eq. (§) yields 



v(.)~^-fi(£)'l>e 



cos 2 7 a-1), (13) 



where cos ja = *a ■ n = sin sin ^ cos (0 — (3 a) + 
cos a. a cos 0. 

The second term in Eq. ( |l3| ) is smaller than the first 
by a factor of order (M/r) 2 — v 4 , and it can be treated 
as a perturbation when integrating the equations of mo- 
tion. Taking the background motion to be a circular or- 
bit in the equatorial plane, we write r(t) = ro[l + eR(t)], 
0(t) = 7r/2 + e6(£), and <j>(t) = Q i + where r and 
O = M 1//2 ro~ 3 / 2 are the background radius and angu- 
lar velocity, respectively, and e = (M/r ) 2 . Substituting 
these relations into the equations of motion derived from 
C and linearizing with respect to e leads to differential 
equations for the unknowns R, O, and Integration of 
these equations is straightforward, and we find 



r(t) = r 



3e 



1+— ^pA^COS 2 UA - 1) 



^^PA sin 2 a a cos2(tt Q t — (3a) 



0(f) = f - T^ ot Xl PA sin 2aA sin (^V - /?/ 



(14) 
(15) 



3c 



1 - — - ^p j4 (3cos 2 a A - 1) 

A 

■ pa sin 2 cos 2(£l i — /3a) 



(16) 



These equations show that the orbits cannot be circular 
in the strict sense, although r and (f> are constant af- 
ter averaging over an orbital period. As we shall see in 
Sec. Ill, the linear growth of 9 signals the occurrence of 
orbital precession. Because we are concerned only with 
the motion over a time scale comparable with the orbital 
period, and because the precessional motion occurs over 
a much longer time scale, we shall ignore this issue here. 

From Eqs. ( |l4] ) and ([16]) we find that the averaged 
orbital radius is given by 



(r) = r 



while the averaged angular velocity is 
3e 



1 - — - ^pa(3cos 2 a A - 1) 



(17) 



(18) 



Eliminating r$ from these equations, we arrive at (we 
discard all terms of order e 2 and higher) 



M 

w 



1/2 



3 ( M 



8 \(r) 



^Pa(3cos 2 a A - 1) 



Inverting this gives 



M 



1 + j v 4 ^pa(3cos 2 ol A - 1) 



(19) 



(20) 



where v = (M(^)) 1 / 3 . 

We shall need also an expression for the orbital energy, 
given by E — \[ix + V(x). Substituting the relations 
(|l4|)-(p0"l) and discarding all terms of order e 2 and higher, 
we find 



E = --(tv* 



1 + ^- v 4 ^2 p A (3 cos 2 a A - 1) 



(21) 



Equations ©-© describe the (averaged) orbital mo- 
tion over time scales comparable to the orbital period. 



III. PRECESSIONS 

It is well known that the Newtonian quadrupole- 
monopole interaction produces a precession of the orbital 
angular momentum, as well as a precession of the spin 
axis of each of the binary companions ^3|. We examine 
this effect here, as well as the precessional motions caused 
by the general relativistic spin-orbit and spin-spin inter- 
actions. 

Calculating the torque associated with the potential 
(|l3| ) and averaging over an orbital period leads to this 
equation describing the orbital precession: 



(L) = n op x L 



(22) 



where the angular velocity of orbital precession is given 
by 



a, 



3 (M 



2 \ r 



Q ~^^PA COSCkA 8 A- (23) 
A 



Here, r = (r) is the orbital radius, and = ((f)) is the 
orbital angular velocity. Equation ( p3| ) shows that the 
time scale for orbital precession is longer than the or- 
bital period by a factor of order (r/M) 2 . Equations ( |2^ ) 
and ( p3| ) also imply that cos aA = L ■ §a does not stay 
constant during the motion. We note that Eq. (|l5|) may 
be recovered by integrating Eq. ( p2|) with L set equal to 
z on the right-hand side. 

On the other hand, treating companion A as a rigid 
body under the influence of the potential (|l^) leads to 
this equation describing the spin precession: 



3 



(s A ) = ttsp.A X 8A, 



(24) 



IV. GRAVITATIONAL WAVES 



where the angular velocity of spin precession is given by 



2 ttia V t ) 



XA 



This shows that the time scale for spin precession is 
longer than the orbital period by a factor of order 
(r/M) 3 / 2 . It should be noted that for simplicity we have 
not accounted for the time variation of cos a a in this 
calculation. 

The quadrupole-monopole interaction is not alone in 
inducing a precession of L and sa- The general relativis- 
tic spin-orbit and spin-spin interactions also contribute 
to these precessions, and this was studied in detail in 
Rcf . [E3 . The spin-orbit contributions to the precessional 
angular velocities are given by 



4mi + 3m2 



and 



2niir 3 



4m2 + 3toi 
2m2r 3 



82 



Ami + 3m2 
2mir 3 



(26) 



(27) 



with a similar equation holding for f2 sP) 2- On the other 
hand, the spin-spin contributions are 



a, 



2r 3 L 

where L = \L\ — /ir 2 f2, and 



(s 2 ■ L) si + (si • L) s 2 



^sp,l — 



1 

2T 3 



s 2 - 3(si • L) L 



(28) 



(29) 



with a similar equation holding for l~2 S p,2- 

From these equations it is easy to see that the spin- 
orbit contribution to the precession dominates over both 
the spin-spin and quadrupole-monopole contributions. In 
fact, 



spin-orbit 



spin-orbit 
spin-spin quadrupole-monopole 



1/2 

T7) • (30) 



The main lesson here is that while the quadrupole- 
monopole precession is entirely a Newtonian effect, it is 
a small one compared with the general relativistic spin- 
orbit precession, and it is of the same order of magnitude 
as the spin-spin precession. 

These precessional motions are important, because 
they produce modulations in the amplitude and phase 
of the gravitational waves emitted by the binary system. 
The modulations produced by the spin-orbit and spin- 
spin interactions were studied in detail in Ref. |^4j . While 
it may be worthwhile to repeat this analysis so as to 
also incorporate the modulations due to the quadrupole- 
monopole interaction, it is doubtful that any of the con- 
clusions would be affected. 



We now use the Einstein quadrupole formula p5[ 



l(Q ab Q ab ) 
5 



(31) 



to calculate the rate of loss of orbital energy to gravita- 
tional radiation, for a binary system moving on a circu- 
lar orbit of (averaged) radius r = (r) and (averaged) 
angular velocity Q = (< t>), where (r) and (4>) are re- 
lated by Eqs. fll9| ) and (P0|). The quadrupole moment 
of the centers of mass is given by Q ab = J^a m A{%A x A ~ 
^rA 2 8 ab ), where, in terms of the relative separation x, 
X\ = — m,2x/M and x 2 = m\x/M. Substituting x = 
r (cos fit, sin fit, 0) into Q ab , and the result into Eq. 
we find 



E 



32 



(32) 



irrespective of the relation between r and f2. Using 
Eq. <M) yields 



5 \M 



1 — v A ^pa(3cos 2 cla — 1) 



(33) 



where v = (MfJ) 1 / 3 . 

The gravitational-wave frequency / is equal to twice 
the orbital frequency, so 



ttM/ = Mfi = v 3 



(34) 



The loss of orbital energy is accompanied by an increase 
in orbital frequency, and therefore, an increase in /. 
The rate of change of the gravitational-wave frequency 
is calculated as / = E/{dE/df) = (3v 2 /irM)E/(dE/dv), 
where E(v) is given by Eq. (pl|). This gives 



• = _96^ u 
J 5nM 3 



1 - -u 4 ^p A (3cos 2 a A - 1) 



(35) 



The coefficient of the 0(v 4 ) term within the square brack- 
ets is the quadrupole-moment contribution to the param- 
eter cr appearing in Eq. (|l|) . Taking into account the def- 
initions (1^) and ([Tl]), this is the same as what was given 
in Eq. (0T 
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